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r£3 Abstract 
Oh 

The physics responsible for neutrino masses and lepton mixing remains unknown. More ex- 
perimental data are needed to constrain and guide possible generalizations of the standard model 
of particle physics, and reveal the mechanism behind nonzero neutrino masses. Here, the physics 
associated with searches for the violation of lepton-flavor conservation in charged-lepton processes 
and the violation of lepton-number conservation in nuclear physics processes is summarized. In the 
^ first part, several aspects of charged-lepton flavor violation are discussed, especially its sensitivity 

to new particles and interactions beyond the standard model of particle physics. The discussion 
concentrates mostly on rare processes involving muons and electrons. In the second part, the sta- 
tus of the conservation of total lepton number is discussed. The discussion here concentrates on 
current and future probes of this apparent law of Nature via searches for neutrinoless double beta 



decay, which is also the most sensitive probe of the potential Majorana nature of neutrinos. 
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1 Introduction 



In the absence of interactions that lead to nonzero neutrino masses, the Standard Model Lagrangian is 
invariant under global U(l) e x U(l)^ x U(1) T rotations of the lepton fields. In other words, if neutrinos 
are massless, individual lepton-flavor numbers - electron-number, muon-number, and tau-number - are 
expected to be conserved. Very concretely, individual lepton-flavor quantum numbers are assigned as 
follows: £ = e, /i, or r, vi have Qi = +1, £ , vi have Qi = —1, while all other fields have charge zero. 
Naturally, the total lepton number associated to U (1) l, the "diagonal" subgroup of U{l) e x U{1)^ x U(1) T 
is a strictly conserved quantity. 

In the last decade a variety of neutrino oscillation experiments proved beyond any doubt that 
neutrino flavors, i.e., individual lepton flavor numbers, are not conserved. After propagating a finite 
distance neutrino beams of a given flavor (electron, muon, tau) contain no longer only neutrinos of the 
initial flavor. That phenomenon has been indisputably established pQ in "disappearance" experiments 
with atmospheric and accelerator neutrinos [y^ -/» and -p> solar neutrinos [y e v e ) and 
reactor neutrinos iv e u e ) and in "appearance" experiments with solar neutrinos iy e — > u^). There 
is also more than three-sigma evidence for the appearance of a new flavor from studies of atmospheric 
(anti)neutrinos (z/ M — > v T ) [2]. 

The only consistent explanation of all of these phenomena is (i) (some of the) neutrino masses are 
nonzero and distinct, (ii) the weakly interacting flavor neutrinos z/ e , z/ M and v T are nontrivial superpo- 
sitions of the so called mass eigenstate neutrinos vi, vi and z/3 (i.e. neutrinos with well-defined masses 
mi, mi and m 3 , respectively). These superpositions are described by a 3 x 3 unitary mixing matrix. 

Further corroborating this picture, very recently, reactor experiments at distances corresponding 
to the "atmospheric oscillation length" [SI IU E] as well as the accelerator "appearance" [y^ — > v e ) 
experiments [HI E] have successfully determined the magnitude of the until now missing mixing angle 
#13. The values of all three mixing angles 612 ~ solar, #23 ~ 6 atm and #13 as well as the two indepen- 
dent mass-squared differences Am^ ~ anc ^ ^ m a*m ~ l^ m iil ~ l^ m 32l are an known now with 
good accuracy (see, for example, [U IH [9] for very recent summaries). Potential CP-violating parame- 
ters in the leptonic mixing matrix, along with the so-called neutrino mass hierarchy, remain virtually 
unconstrained. 

The success of the neutrino oscillation hypothesis is illustrated in Fig. [TJ based on the results of 
the KamLAND reactor oscillation experiment pU]. This textbook- type example clearly shows that 
the initial beam of electron antineutrinos changes into a flavor that is invisible to the detector, i.e., 
oscillates periodically as a function of L/E u , where L is the distance from the neutrino source and E v 
is the neutrino energy. 

The discovery of neutrino oscillations, and its implication of nonvanishing neutrino masses and 
the non-conservation of neutrino flavor, implies that the Standard Model is incomplete and needs 
generalization. The new physics responsible for neutrino masses and mixing remains unknown and will 
only be revealed once more experimental data, from different areas of fundamental physics research, 
become available. Here, we summarize the physics associated with searches for the violation of lepton- 
flavor conservation in charged-lepton processes, and the violation of lepton-number conservation in 
nuclear physics processes. While these two issues seemingly involve different aspects of the physics 
beyond the Standard Model, they are both intimately related to the understanding of the origin of 
neutrino masses, and hence are discussed together. 

The established existence of neutrino flavor violation implies that, barring accidental cancellations, 
all so-called lepton flavor violating processes are also allowed and will occur at some order in perturbation 
theory. The rates for such processes, however, cannot be estimated model-independently and are hence 
expected to provide non-trivial information regarding the nature of new physics. 

On the other hand, neutrino oscillations do not necessarily imply that total lepton number is violated. 
While nonzero neutrino masses and lepton mixing imply that U{l) e x U (1)^ x U (1) T is not a good global 
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Figure 1: Ratio of the measured v e flux to the no-oscillation expectation as the function of 
Lq/E u . Lq = 180 km is the flux weighted average distance from all nuclear reactors. The 
data points are background and geoneutrino subtracted with statistical errors only. The 
curve shows the expectation from neutrino oscillations, taking into account the different 
distances to individual reactors, time variations of the fluxes, and efficiencies. 



symmetry of nature, U(1)l may still remain as an exact global symmetry of the "New Standard Model" 
Lagrangian (one that includes nonzero neutrino masses). This is the case if neutrinos turn out to be 
Dirac fermions. The situation here is analogous to that of the quark sector, where individual quark- 
flavor numbers are known to be violated by weak processes while baryon number U(1)b remains a good 
quantum number (at least as far as the standard model and all current observations are concerned) [j] If 
neutrinos are Majorana fermions, however, U(1)l is not a good symmetry of the New Standard Model. 
Current data are silent concerning the nature - Majorana fermions or Dirac fermions - of the massive 
neutrinos. Different new physics models for nonzero neutrino masses predict the neutrinos to be Dirac 
or Majorana fermions. Only by investigating the validity of lepton-number conservation will we be able 
to distinguish one "type" of new physics from the other. 

The paper is divided into two sections. Section [2] deals with searches for charged-lepton flavor 
violation. Charged-lepton flavor violation (CLFV) is defined as all charged-lepton processes that violate 
lepton-flavor number. These include £ ->■ £'7, £ -> £'£"£"', £ + X ->■ £' + X and X ->■ ££', where 
£,£',... G {e, /1, r} and X are states that carry no lepton-flavor number. We will discuss several aspects 
of CLFV, especially their sensitivity to new particles and interactions beyond the Standard Model. 
We will concentrate on rare processes involving muons and electrons, and will briefly comment on 

1 Both U(1)b and (potentially) U(1)l are only "classical" symmetries of the new Standard Model Lagrangian, being 
violated at the quantum level. U(V)b-l, however, is a potential nonanomalous global symmetry of nature even when 
nonzero neutrino masses are taken into account. 
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processes involving tau-leptons. There are several more comprehensive recent overviews of different 
aspects of CLFV, including [HI [T2J H3J Hi] . Section [3] deals with the total lepton conservation (LNV), 
and in particular with the probe of that law using neutrinoless double beta decay processes. Searches 
for the violation of lepton number often involve nuclear probes and the analyses of the corresponding 
experiments often require deep understanding of nuclear structure. We will briefly review such problems, 
but the main emphasis will be on the possible new physics origins of the violation of this conservation 
law. 

2 Lepton Flavor Conservation 

2.1 Muon to Electron Searches: Current Bounds, Near Future Goals 

Searches for // — > date back to the late 1940's. By the early 1960's it was clear that naive estimates 
for the branching ratio for /i — > cy from the physics responsible for /i — > evv were violated by negative 
experimental results. Indeed, the fact that fi — > ej did not occur with branching ratios above 10 -6 was 
among the most pressing arguments in favor of the so-called 'two-neutrino hypothesis' and the postulate 
that there were two lepton flavors, each accompanied by a conserved lepton-flavor numberj^] In modern 
textbooks, we describe (negative) muon decay as \T — > e'z/gZ/^; the electron neutrino and the muon 
neutrino are different objects. 

Along with /i + — > e + 7 decays, two other rare muon processes provide the most stringent current 
bounds on CLFV, and promise the most near-future sensitivity to these phenomena: fi + — > e + e~e + 
decays, and — > e~ conversion in nuclei. The three are briefly discussed here. Other CLFV processes 
involving muons and electrons include rare kaon decays like Kl — > /i ± e =F , K — > tt// e T , and muonimum- 
antimuonium oscillations, /i + e~ /x~e + pQ. Some of these are very severely constrained and sensitive 
to different types of new physics at different levels. Nonetheless, they will not be considered henceforth. 

fi + — > e + 7 decays are currently being pursued by the MEG experiment at the Paul Scherrer Institut 
([15J and references therein). Presently, Br(/z + — > e +r y) < 2.4 x 10~ 12 at the 90% confidence level [To] . 
MEG aims at being sensitive to Br(yU + — > e +r y) > 1CT 13 , while potential improvements may lead to 
sensitivities approaching 10~ 14 [16jj^] Experimentally, it is extremely challenging to envision improving 
the sensitivity to Br(/x — > ej) beyond 10~ 14 (for a discussion, see [12]) - the current sensitivity is already 
dominated by accidental backgrounds which grow with the muon beam intensity. For this reason, there 
are no proposals for improving on the sensitivy to ji — > beyond what can be ultimately achieved by 
MEG. 

Presently, the Br(/i + — > e + e~e + ) < 1.0 x 10~ 12 at the 90% confidence level [18]. This limit was 
obtained over 20 years ago. Presently, a letter of intent has been submitted to the Paul Scherrer 
Institut to conduct an improved search for fi + — > e + e~e + , ultimately sensitive to branching ratios 
above 10~ 15 or, perhaps, 10~ 16 [19]. As with /i + — > e + 7 decays, sensitivity beyond 10~ 16 , regardless of 
the intensity of the muon beam, appears very challenging. 

— > e~ conversion in nuclei is the process where a bound muon interacts with the nucleus and 
converts into an electron, which is born with enough kinetic energy to "escape" the Coulomb potential: 
Ijl~N — >■ e~N ( -*\ where N is some nucleus. Experimentally, the signal for — > e~ conversion in nuclei is 
a monochromatic electron, whose energy lies just beyond the kinematical end-point of Michel electrons 
produced by muon decay in orbit. The — > e~ conversion rate is usually expressed in units of the 
capture rate, fi~N — > u^N'. Presently, the most stringent bound is Br(/i — > e conv in Au) < 7 x 10~ 13 
at the 90% confidence level [20] ■ There are a few proposals aimed at experiments sensitive to Br(/i — >■ 

2 The hypothesis was expanded in the 1970's to include a third conserved lepton-flavor number. 

3 After the completion of this review, the MEG experiments made public their most recent result, Br(/i + — > e + 7) < 
5.7 x 1CT 13 at the 90% confidence level [17]. 
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e conv m Al) > lCT 16 or betteiQby the end of this decade j2H I22J- Unlike the two rare decays discussed 
earlier, searches for — > e~ conversion in nuclei are not, naively, expected to hit any experimental 
"wall" until conversion rates below 1CT 18 or lower [12]. Experimental setups sensitive to conversion rates 
below 10~ 17 are currently under serious study both at Fermilab (assuming Project X becomes available) 
and J-PARC. Experimentally, in the long-run, it is widely anticipated that fi~ — > e~ conversion in nuclei 
will provide the ultimate sensitivity to CLFV. 



2.2 v Standard Model Expectations 

In spite of the fact that we have determined that CLFV must occur, measurements of neutrino oscillation 
processes do not allow us to reliably estimate the rate for the various CLFV processes. The reason is 
that while neutrino oscillation phenomena depend only on neutrino masses and lepton mixing angles, 
the rates for the various CLFV processes depend dramatically on the mechanism behind neutrino masses 
and lepton mixing, currently unknown (for recent reviews see, for example, [23]). Different neutrino 
mass-generating Lagrangians lead to very different rates for CLFV. Some of these will be discussed 



briefly here and in Sec. 2.4 



The massive neutrino contribution to CLFV that involves only active neutrinos is absurdly small. 
For example [23] . 

Ml 



Br(/i ->• erf) 

62,71 



i=2,3 w 



<1(T 54 , (1) 



where U a i are elements of the neutrino mixing matrix, Am^ are the neutrino mass-squared differences, 
a is the fine-structure constant, and Mw is the VT-boson mass. Similar ridiculously small rates are 
expected for /i — > eee, // — > e conversion and rare process involving taus. The estimate above applies 
to some neutrino mass models, including minimal scenarios with Dirac neutrinos. The reason behind 
the tiny branching ratios is well-known. CLFV, as defined above, is a flavor- changing neutral current 
process and such processes are subject to the GIM mechanism (or generalizations thereof). 

In many neutrino-mass generating scenarios, the active neutrino contribution turns out to be, not 
surprisingly, severely subdominant. In the famous seesaw mechanism [251 [26] . for example, heavy 
neutrino contributions to CLFV are naively expected to be of order those of the light neutrinos, but there 
is no theorem that prevents them from being much, much larger. According to [27] , for example, current 
experimental constraints on the seesaw Lagrangian allow Br(r — > /ry) as large as 10 -9 , Br(/i — > ej) 
as large are 4 x 10~ 13 , and normalized rates for ft — > e conversion in nuclei that saturate the current 
experimental upper bound. For more details, see, for examples, [28J. 



2.3 Some Model Independent Considerations 

Independent from the mechanism behind neutrino masses, it is often speculated that the rates for 
different CLFV processes are, perhaps, just beneath current experimental upper bounds. The reason 
is we suspect, for several reasons, that there are new degrees of freedom beyond those in the Standard 
Model. We also suspect that some of those have masses around 1 TeV. Since lepton-flavor numbers are 
known not to be good quantum numbers, it is generically expected that virtual processes involving the 
new degrees of freedom will mediate, at some order in perturbation theory, CLFV. 

Some concrete new physics scenarios will be briefly discussed in Sec. [Z4 Here we discuss two 



effective Lagrangians^Jthat mediate CLFV processes involving muons aiming at illustrating how searches 
for CLFV are sensitive to new physics, and how different CLFV channels compare with one another. 



4 Some proposals involve other nuclei, including 48 Ti. 

5 Parts of this discussion were first presented in writing in pilfTB] . 
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After integrating out heavy degrees of freedom, and after electroweak symmetry breaking, CLFV 
is mediated by effective operators of dimension five and higher. We first concentrate on the following 
effective Lagrangiarj^] 

TTl 

£clfv = ^ma ^ 1 ^ + hx ' 

Hl^L {u L Yu L + d L ^d L ) + h.c. . (2) 



/V 



(1 + k)A 2 

The subscripts L, R indicate the chirality of the different Standard Model fermion fields, is the 
photon field strength and m M is the muon mass. The coefficients of the two types of operators are 
parameterized by two independent constants: the dimensionful A parameter (with dimensions of mass), 
which is meant to represent the effective mass scale of the new degrees of freedom, and the dimensionless 
parameter k, which governs the relative size of the two different types of operators. The magnetic- 
moment type operator in the first line of Eq. ^ directly mediates 11 — > e'y and mediates /i — > eee and 
fi — > e conversion in nuclei at order a. The four-fermion operators in the second line of Eq. (|2]), on the 
other hand, mediate \i — > e conversion at the leading order and [i — > ej, ii — > eee at the one-loop level. 
For k <C 1, the dipole-type operator dominates CLFV phenomena, while for k ^> 1 the four-fermion 
operators are dominant. 

The sensitivity to A as a function of k for ii — > ej and \i — > e conversion efforts is depicted in Fig. [2j 
For k < 1, an experiment sensitive to Br(/i — > cy) > 10~ 13 will probe A values less than 2500 TeV, 
while for k > 1 an experiment sensitive to Br(/x — > e conv in 27 A1) > 10~ 16 will probe A values less 
than 7000 TeV. 

Relevant information can be extracted from Fig. [2j CLFV already probes A values close to 1000 TeV 
and next-generation experiments will start to probe A ~ 10 4 TeV and beyond. Furthermore, a /i — > 
e conversion experiment is "guaranteed" to outperform a /i — > cy experiment for any value of k as 
long as it is a couple of orders of magnitude more sensitive. Since, as already discussed, it appears 
very challenging to perform a fi — > ej experiment sensitive to branching ratios smaller than 10~ 14 , 
fi — > e conversion searches (not expected to hit any "wall" before normalized rates around at most 
10~ 18 ), are the more effective way of pursuing CLFV after the on-going MEG experiment is done 
analyzing its data. 

Similarly, we can ask what are the consequences for CLFV if the new physics is best captured by 
the following "leptons-only" effective Lagrangian: 

TTl 

£clfv = (ft. + ^^ ^rv^lF^ + h.c. 

^Ll^e L (e7 M e) + h.c. . (3) 



(1 + k)A 2 

Similar to the dimension-six operators in the second line of Eq. ([2]), the dimension-six operator in the 
second line of Eq. ^ mediates ii — > eee at the tree level and ii — > e^y, fx — > e conversion at the one- 
loop level. Similar to Eq. ([2]), the dimensionless parameter k determines whether the dipole-like or the 
four-fermion interaction is dominant when it comes to CLFV. 

The sensitivity to A as a function of k for ii — > e'j and ii — > eee efforts is depicted in Fig. [3] Here, 
for k ^> 1, an experiment sensitive to Bi(fi — > eee) > 10~ 15 will probe A values less than 1800 TeV. As 
in the example depicted in Fig. [2j we note that a ii — > eee experiment is guaranteed to outperform a 
ix — \ e7 experiment, for any value of k, as long as it is a few hundred times more sensitive. Whether 
this can be realistically achieved in future experiments is still under investigation [121 HH]- 



3 The most general effective Lagragian includes several other terms [14. The subsets included in Eqs. (2j3), however 



are sufficient to illustrate all issues discussed here. Modulo extreme constructive/destructive interference effects among 
different effective operators, the points made here remain valid. 
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Figure 2: Sensitivity of a // — > e conversion in 27 A1 experiment that can probe a normalized 
capture rate of 10~ 16 and 10~ 18 , and of a fi — > e^ search that is sensitive to a branching ratio 
of 10 -13 and 10 -14 , to the new physics scale A as a function of k, as defined in Eq. (pi). Also 
depicted is the currently excluded region of this parameter space. 



A model independent comparison between the reach of /i — > eee and [l — > e conversion in nuclei is 
a lot less straight forward. If the new physics is such that the dipole-type operator is dominant (k < 1 
in Figures [2] and [3]), it is easy to see that near-future prospects for [i — >■ e conversion searches are 
comparable to those for fi — > eee, assuming both can reach the 10 -16 level, fi — > e conversion searches 
will ultimately dominate, assuming these can reach beyond 10 -17 , and assuming /i — > eee searches 
"saturate" at the 10 -16 level. Under all other theoretical circumstances, keeping in mind that k and A 
are not the same, it is impossible to unambiguously compare the two CLFV probes. 

The discussions above also serve to illustrate another "feature" of searches for CLFV violation. 
In the case of a positive signal, the amount of information regarding the new physics is limited. For 
example, a positive signal in a /i — > e conversion experiment does not allow one to measure either A or 
k but only a function of the two. In order to learn more about the new physics, one needs to combine 
information involving the rate of a particular CLFV process with other observables. These include other 
CLFV observables (e.g., a positive signal in /i — > ej and /i — > eee would allow one to measure both 
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Figure 3: Sensitivity of a fi — > eee experiment that is sensitive to branching ratios 10~ 14 
and 10~ 16 , and of a // — > search that is sensitive to a branching ratio of 10 -13 and 10 -14 , 
to the new physics scale A function of clS defined in Eq. Also depicted is the 
currently excluded region of this parameter space. 

k and A if Eq. ^ describes CLFV), studies of electromagnetic properties of charged leptons (g — 2, 
electric dipole moments), precision studies of neutrino processes (including oscillations), and, of course, 
"direct" searches for new, heavy degrees of freedom (Tevatron, LHC). Valuable information, including 
the nature and chirality of the effective operators that mediate CLFV, can be obtained by observing 
\i — y e conversion in different nuclei [El [29], [30] or by studying the kinematical distribution of the 
final-state electrons in fi — > eee (see [H] and references therein). 

Before moving on to specific new physics scenarios, it is illustrative to compare, as model-independently 
as possible, new physics that mediates CLFV and the new physics that may have manifested itself in 
precision measurements of the muon anomalous magnetic moment. In a nutshell, the world's most 
precise measurement of the g — 2 of the muon disagrees with the world's best Standard Model estimate 
for this observable at the 3.6a level (for an updated overview see pQ, and references therein). New, 
heavy physics contributions to the muon g — 2 are captured by the following effective Lagrangian: 

C g -2 D -rf /^<V/i^ + h.c. . (4) 
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Current g — 2 data point to A ~ 8 TeV. Eq. Q, however, is very similar to Eqs. (2j3) in the limit 
k < 1, keeping in mind that A in Eq. Q need not represent the same quantity as A in Eqs. (2]3) in 
the limit k< 1. 

We can relate A in Eq. Q to that in Eqs. (2]3) in the following suggestive way: (Aclfv) 2 = 
9 eil (A g -2) 2 - Here the parameter 6 etl measures how flavor-conserving is the new physics. For example, 
if # eAt = 0, the new physics is strictly flavor conserving, while if the new physics is flavor- indiferent, 
efl ~ 1. In the latter case, negative searches for \i — > already preclude a new physics interpretation 
to the muon g — 2 results as these constrain A > 1000 TeV. On the other hand, if the muon g — 2 
discrepancy is really evidence for new physics, searches for fi — > ej reveal that the "amount" of flavor 
violation in the new physics sector is very small: 8 efJ- < 10~ 4 . A more detailed comparison of these two 
probes of new physics can be found, for example, in [31J. 



2.4 Some Scenarios and Examples 

Models for new physics that introduce new degrees of freedom at the weak scale are often severely 
constrained by, or make specific predictions for, CLFV observables. These include new physics scenarios 
that address the gauge hierarchy problem and certain models for nonzero neutrino masses. Concrete 
predictions further require that the models specify the flavor structure of the new physics, and are also 
obstructed by the fact that some of the flavor parameters of the lepton sector, including the neutrino 
mass hierarchy and the value of all CP-violating parameters, remain unknown. Here we discuss, in a 
little detail, a couple of examples aimed at conveying the most important lessons and challenges, and 
list other possibilities. 

Neutrino Majorana masses can be generated if the Standard Model particle content is augmented 
to include a Higgs boson SU(2) triplet T with hypercharge one [26]. In this case, the gauge symmetries 
allow 

C D XapLaTLp + h.c, (5) 

where \ a p = \p a are dimensionless Yukawa couplings, L a are the SU{2) lepton doublet fields, a,(3 = 
e, (jl, t. After the neutral component of T acquires a vacuum expectation values vt, neutrinos acquire a 
Majorana mass matrix m u = Xvt- 

All CLFV processes are mediated by one-loop diagrams involving the Higgs triplet fields, while 
£ — » £'£"£"' processes are mediated at the tree level via the exchange of the doubly-charged components 
of T (for details see, for example, [321 EE] and references therein). This scenario manifests itself at 
low-energies in the form of (a variant of) Eq. ^ with, in general, large k. 

The branching ratio for fi — > eee, for example, is 
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>-_ r — ( 6) 

where is the mass of the doubly-charged component of T. The rates for other CLFV observables 
are also functions of vt, Mt, and the elements of m v . Concrete predictions for the rates of CLFV 
processes, in this case, are not possible without more external assumptions; there is no input regarding 
Mt and vt (unless triplet Higgs bosons are discovered at the LHC!). Correlations among different CLFV 
observables, however, are well-defined; m u is, in principle, uniquely determined by neutrino masses and 
lepton mixing observables. Recent detailed studies, which take into account the recent discovery that 
#13 is not zero, can be found in (28] [33] . 

A class of very well studied models that usually predicts k 1 is that of supersymmetric versions 
of the Standard Model with R-parity conservation. Superpartner loops involving gauginos and sleptons 
lead to [33], in the language of Eqs. ( 2p ), 

— g2e fl m 

A 2 ~ 16tt 2 M| usy e/ " 1 ' 
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where e, g are the electromagnetic and SU(2) L couplings, respectively, and M SU sy is some effective 
superpartner mass. The flavor violating parameter 9 efl can be expressed in terms of the off-diagonal 



,2 



slepton mass-squared matrix element Am^. 

Am 2 







lie 

JW SUSY 



Similar to the g — 2 discussion in the previous section, weak-scale superpartner masses are only allowed 
by CLFV data if 6* Me is small, which, in this case, translates into Am~g much smaller than the diagonal 
slepton masses-squared. 

One way to guarantee that this is the case is to assume that the mechanism of supersymmetry 
breaking, which determines the slepton mass-squared matrix, is flavor blind so that, in the absence 
of Standard Model flavor effects, Am~ g = 0. Since, in the Standard Model with nonzero neutrino 
masses, flavor-numbers are not conserved, quantum mechanical effects lead to nonzero, calculable Arn-g. 
Concrete predictions, alas, depend dramatically on the mechanism behind neutrino masses. In the case 
of the Type-I seesaw mechanism, for example, 

k k 

where y are the neutrino Yukawa couplings, M^ k are the Majorana masses of the k — 1, 2, 3 right-handed 
neutrinos and Muv is the ultraviolet scale where the supersymmetry breaking parameters originate]^] 

Predictions for CLFV processes will depend, in this case, on the values of the neutrino Yukawa 
couplings, which are not directly constrained by neutrino mass and lepton mixing data. If these are large, 
expectations are that the rates for several CLFV processes, especially // — > e^j and fi — > e conversion 
in nuclei, are large enough to be seem in current or next-generation experiments. For recent detailed 
study which discusses different ansatze for the neutrino Yukawa couplings see, for example, [351 [37]. 

Equation ^ reveals that, at least in principle, one can relate the rates for CLFV processes to a par- 
ticular combination of neutrino Yukawa couplings. On the other hand, neutrino masses reveal, through 
the seesaw mechanism, a different combination of neutrino Yukawa couplings: m u oc Y^kV^M^ypk- 
The hope is that by combining neutrino masses, CLFV and the measurement of super-partner masses 
and couplings at high energy colliders one can determine the neutrino Yukawa couplings and the right- 
handed neutrino masses. This would not only nail down the physics behind neutrino masses but would 
also allow one to test whether the baryon asymmetry of the universe arises due to leptogenesis (see, for 
example, [38] and references therein). Needless to say, in order to test high scale leptogenesis, several 
"minor miracles" have to occur, including: there must be weak scale supersymmetry, and we must 
measure several supersymmetry breaking parameters very precisely; we must understand the mecha- 
nism of supersymmetry breaking reliably; we need to observe several CLFV events in several channels 
(including some involving tau-leptons); and a few others. For a detailed discussion see, for example, 
[39]. 

Other new physics scenarios will lead to different values of k, including k ~ 1 and k ^> 1. Supersym- 
metric models with trilinear R-parity violation serve as an excellent laboratory for exploring different 
possibilities (30] E] . In the presence of, say, large couplings A' for the LQD C superpotential term (see 
[30] for details and notation), one realizes Eq. ^ in the limit k>1 and 

— ~ -± (10) 

A 2 M 2 ' 1 ; 

In such a scenario the rate for \i — > e conversion in nuclei is much larger (several orders of magnitude) 
than that for — y e^y and fi — > eee. The presence of other R-parity violating couplings may lead to 



7 A different example, where neutrino masses originate from the Type-II seesaw, is discussed in detail 
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scenarios where k ~ 1 both in Eq. (J2J) and Eq. (J3J), while yet another choice leads predominantly to 
Eq. ^ in the limit K > 1. 

Large rates for CLFV are not a privilege of supersymmetric scenarios. Any complete new physics 
scenario that addresses the gauge hierarchy problem will lead to CLFV rates that are much, much larger 
than those associated to the naive massive neutrino contribution (see Eq. ([!])). Concrete examples, 
associated to different values of k, have been calculated for models with large extra dimensions [42J, 
Randall-Sundrum models with fermion fields in the bulk [13], little Higgs models jH], and many others. 

2.5 CLFV with Tau Leptons 

The most stringent bounds on tau-lepton-flavor violation come from rare tau decay processes. Thanks 
to the B-factories, branching ratios for r — > £q, r — > £££' and r — > £ + X, where £, £' = e, \i and X are 
different quark states, are known to be, for the most sensitive decay modes, below a few times 10 -8 pp. 
Super-B factories, which aim at integrated luminosities larger than 10 ab^ 1 , aim at being sensitive to 
branching ratios larger than 10~ 9 (or even a little lower), especially when it comes to "background-free" 
decay modes (e.g. r — > e + /i _ 

It is not possible to compare, in a model independent way, the sensitivity of rare tau and rare muon 
processes. At face value, rare muon processes are at least a few orders of magnitude more sensitive than 
rare tau processes. This remains true when one compares qualitatively similar processes, e.g., // — > eq 
versus r — > eq, fi — > eee versus r — > efifi, 11 — > e conversion in nuclei versus r — >• ep°, etc. 

There are several new physics scenarios that predict significantly enhanced rare tau process compared 
to rare muon processes. Enhancement factors depend, of course, on the model. These include chirality 
effects, which tend to enhance the rate of tau processes by a factor (m T /m M ) 2 ~ 280 and different effects 
related to neutrino masses or the pattern of lepton mixing. For example, r — > [iq may be enhanced 
with respect to fi — > eq by a factor Am^ 3 /Amf 2 ~ 30 or \U r3 /U e3 \ 2 ~ 20. 

2.6 Final Thoughts 

The discovery of neutrino oscillations reveals that charged-lepton flavor violating phenomena must 
occur. Naive massive-neutrino expectations are that the rates for CLFV processes are absurdly small 
thanks to the GIM mechanism and the fact that neutrino masses are tiny compared to the weak scale. 
Nonetheless, if there are new degrees of freedom at the TeV scale, rates for CLFV processes are expected 
to be much larger. In most complete scenarios that address the gauge hierarchy problem by introducing 
new degrees of freedom around the weak scale, expected rates for CLFV are very large, and one is often 
required to add new conditions or constraints to the new physics sector in order to satisfy bounds from 
negative searches for CLFV. 

In some new physics scenarios, fi — > eq is the most promising channel for pursuing searches for 
CLFV. This, however, is not a universal statement. Indeed, searches for fi — > e conversion prove to 
be the ultimate probes of CLFV. Regardless, all searches for CLFV in all channels should be pursued. 
This is especially true if CLFV is observed in some channel. By combining the results from different 
probes we will be able to learn much more about the nature of the new physics. 

Finally, CLFV phenomena may be intimately related to the physics behind neutrino masses. They 
may also play a key role in our understanding of the seesaw mechanism, grand unified theories, and 
the physics behind the matter-antimatter asymmetry of the universe. Regardless of what new (or old) 
physics we may discover at the LHC, future studies of CLFV phenomena are guaranteed to provide 
unique, complementary information - even if the LHC fails to detect any new degrees of freedom at the 
TeV scale. 
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3 Lepton number conservation 



3.1 Origin of the neutrino mass 

Neutrino masses are six or more orders of magnitude smaller than the masses of other charged fermions. 
Moreover, the pattern of masses, i.e. the mass ratios of neutrinos, is rather different (even though it 
remains largely unknown, one does not know whether it follows the normal or inverted hierarchy or 
whether the pattern of masses is essentially degenerate ) than the pattern of masses of the up- or down- 
type quarks or charged leptons. And the neutrino mixing matrix does not look at all as its analog, the 
CKM mixing matrix of the down- type quarks. All of these facts suggest that, perhaps, the origin of 
the neutrino mass is different than the origin (which is still not well understood but is believed to be 
the coupling to the Higgs vacuum expectation value) of the masses of the other charged fermions. 
The smallness of the neutrino masses can be possibly understood following the finding of Weinberg 
who pointed out more than thirty year ago that there exists only one lowest order (dimension 5, 
suppressed by only one inverse power of the corresponding high energy scale A) gauge-invariant operator 
given the content of the Standard Model. After the spontaneous symmetry breaking when the Higgs 
acquires vacuum expectation value that operator represents the neutrino Majorana mass which violates 
the total lepton number conservation law by two units, 

C™ = ^(i?v)+h.c. , (11) 

where v ~ 250 GeV and might be of the order of unity. For sufficiently large scale A neutrino 
masses are arbitrarily small. There exist a large number of physics realization of this basic idea. They 
all necessarily lead to the conclusion that neutrinos are Majorana fermions and that the total lepton 
number is not exactly conserved quantity. 

The most popular physics model explanation of the smallness of neutrino mass is the seesaw mecha- 
nism, which is also roughly thirty years old [251 126]. In it, the existence of heavy right-handed neutrinos 
Nr is postulated, and by diagonalizing the corresponding mass matrix one arrives at the formula 

= 1 (12) 

where the Dirac mass mn is assumed to be a typical charged fermion mass and is the Majorana 
mass of the heavy neutrinos Nr. Again, the small mass of the standard neutrino is related to the large 
mass of the heavy right-handed partner. Requiring that m u is of the order of 0.1 eV means that 
(or A) is ~ 10 14-15 GeV, i.e. near the GUT scale. That makes this, so-called Type I seesaw template 
scenario particularly attractive. However, clearly an observation of the gauge-singlet neutrinos Nr, if 
they exist, is not possible. 

There are, however, other possible scenarios that lead to a small neutrino Majorana mass. In the 
Type II seesaw one adds a Higgs triplet (£ ++ ,£ + ,£°) which couples directly to the symmetric triplet 
combination of two (/, v)l doublets. In that case 

m„ = M£°> , (13) 

that could be small provided that the is very small, even for a "natural" value of the Yukawa 
coupling h v (see the discussion following Eq. (|5]) in Section IID). It works since the spontaneous 
breaking of electroweak symmetry is accomplished by the vacuum expectation value of (<&°) = v , of the 
ordinary Higgs and thus could be small. If the mass of the Mg ~ 1 TeV then its decay into 
l^lj could be observable and would be another signal of the lepton number violation. In Type III the 
singlet neutrinos Nr are replaced by a fermion triplet S + , S°, X - . Small neutrino Majorana masses are 
again obtained. (See Ref. [16] for a pedagogical treatment of the seesaw mechanism.) 
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Figure 4: By adding loops involving only standard weak interaction processes the Ou/3/3 
decay amplitude (the black box) implies the existence of the Majorana neutrino mass [4"T] . 



All of these possibilities not only offer an explanation of the smallness of neutrino masses, but 
they also suggest that neutrinos are likely massive Majorana fermions. To test such hypothesis one 
should explore various ways of lepton number violation. In that context it is important to point out, 
as shown by Schechter and Valle [47] . that the observation of lepton number violating process, in the 
case considered by them, the neutrinoless double beta decay dd — > uuee, implies that neutrinos are 
massive Majorana particles (see Fig. [4] for illustration). This theorem, however, does not mean that 
it is easy (or even generally possible) to determine the neutrino mass once the 0^/3/3 decay is observed 
and its decay rate determined. In fact, the multiloop graph in Fig. [4] only illustrates the fundamental 
relation between the Oz//3/3 decay and the Majorana neutrino mass. Taking into account the upper limit 
for the Ov(3(3 amplitude, deduced from the lower limit of the Oz//3/3 decay half-life, one comes to the 
conclusion that the Feynman graph in Fig. [4] corresponds to an unobservably small neutrino Majorana 
mass jlH]. The true origin of the neutrino masses must be elsewhere, probably in one of the possible tree 
level diagrams. Nevertheless the graph in Fig. [4] illustrates the intimate relation between the neutrino 
Majorana mass on one hand, and the neutrinoless double beta decay on the other hand. 



3.2 Testing the total lepton number 

As stresses in the Introduction, testing whether the total lepton number is conserved, i.e. whether 
U(1)l is a good symmetry or not, is of fundamental importance. It separates the form of the "New 
Standard Model" into two distinct classes with massive Dirac or Majorana neutrinos. 

There are various ways to test whether the total lepton number conservation. No violation has been 
seen so far. Examples of the potentially lepton number violating (LNV) processes with important limits 
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are 



(Z, A) -> (Z + 2, A) + 2e~; half - life > 10 25 years 

yT + (Z, A) — > e + + [Z — 2, A); exp. branching ratio < 10 

K + —7- / u + /i + 7r _ ; exp. branching ratio < 3 x 10~ 9 , 

v e emission from the Sun; exp. branching ratio < 1CT 4 . 
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(14) 



Detailed analysis suggests that the study of the 0v/3/3 decay, the first on the list above, is by far 
the most sensitive test of LNV. In simple terms this is caused by the amount of tries one can make. A 
100 kg Oz//3/5 decay source contains ~ 10 27 nuclei that can be observed for a long time (several years). 
This can be contrasted with the possibilities of first producing muons or kaons, and then searching for 
the unusual decay channels. The Fermilab accelerators, for example, produce ~ 10 20 protons on target 
per year in their beams and thus correspondingly smaller numbers of muons or kaons. However, the 
statement about the overwhelming sensitivity to LNV of the Ovflfl decay is based on the assumption 
that the source of the lepton number nonconservation is the fact that the light neutrino mass eigenstates 
ui, V2 and ^3 are selfconjugate Majorana particles. It also assumes that there is no symmetry that would 
cause the corresponding m ee entry of the neutrino mass matrix in the flavor basis to vanish. 

However, the origin of the lepton number nonconservation could be also elsewhere. It can be related 
to the existence of ~ TeV mass, so far unobserved particles, e.g. a double charged Higgs £ ++ of the 
seesaw Type II scenario, that can directly decay into a pair of the same charge leptons. If one could 
produce a sufficient number of such particles, and detect their decay, that would represent another test 
of the lepton number conservation. Examples of the particle physics models of such type and estimates 
of the corresponding production cross sections could be found e.g. in Ref.|49j where other references 
for analogous ideas could be found. Some of these processes could, in principle, lead to the observation 
of LNV even given the existing limits from the Oz//3/3 decay. 

Nevertheless, in the following we shall discuss the search for LNV using the 0u(3(3 decay as the probe, 
and the corresponding physics involved. 

3.3 Brief history of 0v(3(3 decay searches 

In double beta decay two neutrons, bound in the ground state of an even-even initial (or parent) 
nucleus are transformed into two bound protons, typically also in the the ground state of the final (or 
granddaughter) even-even nucleus, with the simultaneous emission of two electrons only for the Ov/3/3 
mode, or two electrons plus two v e for the 2u(3(3 mode 



The (3(3 decay, in either mode, can proceed only if the initial nucleus is stable against the standard 
(3 decay (both (3~ and (3 + or EC). That happens exclusively in even-even nuclei where, moreover, 
the ground state is always I n = + . The (3(3 decay rate is a steep function of the energy carried by 
the outgoing leptons (i.e. of the decay Q- value). Hence, transitions with larger Q-value are easier 
to observe and experimental search is centered on such nuclei. Possible exception to this rule is the 
resonant neutrino less double electron capture (OuECEC) where Q — > 0. This topic will not be further 
discussed here, but its theory is comprehensively covered in, e.g. Ref. [50] . 

In both modes of the (3(3 decay the rate can be expressed as a product of independent factors that 
depend on the atomic physics (the so called phase-space factors G 0u and G 2u ) that include also the 
Q-value dependence as well as the fundamental physics constants which can be evaluated accurately, 
nuclear structure (the nuclear matrix elements M 0u and M 2p ) that need to be evaluated using the 
nuclear structure theory, and thus are known but only with considerable uncertainty, and for the Ou(3(3 



(Z, A) g . B . ^(Z + 2, A) g _ s . + 2e- + (2z/ e ) . 



(15) 
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Figure 5: Moore's law of 0^/3/3 searches. The past experiments are characterized by the 
sensitivity to the effective neutrino Majorana mass (mpp). 



mode the possible particle physics parameters (the effective neutrino mass (rripp) in the simplest case). 
These particle physics parameters are, naturally, the important output of the search for the 0uj3j3 decay. 
Thus 

1 = c^M^Km^ ; = G^|M 2 f . (16) 

J l/2 J l/2 

The phase space factors G 0u and G 2u can be found in More recent, somewhat different values 
of them that include effects of the nuclear size, are published in [52J. 

The search for the Ozv/3/5 decay has a very long history. There have been over hundred published 
results in the last fifty years. The most sensitive ones, expressed as the upper limits of the effective 
Majorana mass < mpp > are shown in Fig. [5j Over most of that time the trend of the time dependence 
more or less followed a straight line on the semilog plot as indicated by the blue dashed line (the 
Moore's law of (3/3 decay). That led to the improvement in sensitivity to < rripp > by about an order of 
magnitude per decade (i.e. by a factor of ~100 in half-life). However, it is a bit worrisome to note that 
during the last decade that trend was no longer followed, and only slight improvement was achieved 
very recently. This shows that the experiments have become qualitatively more complex and costly; it 
might be challenging to maintain the historic trend in near future. 

At the same time, during the last two decades the less fundamental 2v(3(3 decay has been observed 
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in "live" laboratory experiments in many nuclei, often by different groups and using different methods. 
That shows not only the ingenuity of the experimentalists who were able to overcome the background 
nemesis, but makes it possible at the same time to extract the corresponding 2v nuclear matrix element 
from the measured decay rate. Study of the 2v(3(5 decay is an interesting nuclear physics problem 
by itself. Moreover, evaluation of the M 2u matrix elements is an important test for the nuclear theory 
models that aim at the determination of the analogous but different quantities for the more fundamental 
Oz/ neutrinoless mode. (We will discuss that issue later.) The resulting nuclear matrix elements M 2u , 
which have the dimension energy -1 , are plotted in Fig. [6] (The corresponding experimental half-lives 
Ti /2 were taken from the Refs. [53l EH [55l [561 EZl EH! ESI [601 ED E2].) Note the pronounced shell 
dependence; the matrix element for 100 Mo is almost ten times larger than the ones for 130 Te or 136 Xe. 



0.25 




0.15 



0.05 



48 76 82 96 100 116 130^ 136 v 150 XTJ 

Ca Ge Se Mo Cd Te Xe Nd 



Figure 6: Matrix elements M 2v in MeV 1 based on the experimental halflife measurements. 



3.4 Mechanism of the Ov/3/3 decay 

The relation between the O^/3/3-decay rate and the effective Majorana mass (mpp) is to some extent 
problematic. The rather conservative assumption leading to Eq.(16) is that, if Qv(3/3 decay will occur 



at an observable rate, it will proceed dominantly through the exchange of a virtual light, but massive, 
Majorana neutrino between the two nucleons undergoing the transition, and that these neutrinos interact 
by the standard left-handed weak currents. But that is not the only theoretically possible mechanism. 
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Lepton number violating (LNV) interactions involving so far unobserved much heavier (~ TeV) particles 
might potentially lead to a comparable 0uj3j3 decay rate. 

In general Ou/3/3 decay can be generated by (i) light massive Majorana neutrino exchange, as in the 
Eq.(16) or (ii) heavy particle exchange (see, e.g. Refs.[63| IM]). resulting from LNV dynamics at some 
scale A above the electroweak one. (That can happen e.g. in the seesaw Type II and III scenarios.) The 
relative size of heavy (Ah) versus light particle (A£) exchange contributions to the decay amplitude 
can be crudely estimated as follows [6"5] : 

A L ~Gl^- A H ~ Gi— ^ ~ M ^ k2) (17) 
L (k 2 ) ' H A 5 ' A L A 5 (m^} ' 1 ' 

where (mpp) is the effective neutrino Majorana mass, (k 2 ) ~ (100 MeV) 2 is the typical light neutrino 
virtuality, and A is the heavy scale relevant to the LNV dynamics. Therefore, Ah/Al ~ 0(1) for 
(mpp) ~ 0.1 — 0.5 eV and A ~ 1 TeV, and thus the LNV dynamics at the TeV scale would lead to similar 
Oi//3/3-decay rate as the exchange of light Majorana neutrinos with the effective mass (tripp) ~ 0.1 — 0.5 
eV. 

Obviously, the 0v(3/3 lifetime measurement by itself does not provide the means for determining the 
underlying mechanism. The spin-flip and non-flip exchange can be, at least in principle, distinguished by 
the, naturally very difficult, measurements of the single-electron spectra or of the electron polarizations. 
However, in most cases the mechanism of light Majorana neutrino exchange, and of heavy particle 
exchange, cannot be separated by the observation of the emitted electrons. 

The main difference between the mechanism involving the light massive Majorana neutrino exchange 
and the heavy (~ TeV) particle exchange is the range of the operator that causes the transition. The 
light neutrino exchange represents two point-like vertices separated by the distance r ~ 1/q, with 
q ~ O(100) MeV. The decay rate is then proportional to the square of the effective Majorana neutrino 
mass as in Eq. ( 16 ). On the other hand the heavy particle exchange represents a single point-like vertex 
(six fermions, four hadrons and two leptons), i.e. dimension 9 operator. Proper treatment of the short 
range nucleon-nucleon repulsion and the effects related to the finite size of the nucleon is obviously 
crucial in that case. The relation between the neutrino mass and the decay rate is very indirect in that 
case, but extraction of various important particle physics parameters is possible when the corresponding 
particle physics model has been specified. 



3.5 Specific models: LRSM and RPV-SUSY 

Lets consider the case where the scale of LNV is relatively low (~ TeV). The relevant particle physics 
models in that case usually contain not only the possibility of Lepton Number Violation but also a 
possibility of the Lepton Flavor Violation (CLFV) involving charged leptons. Denoting the new physics 
scale by A, one has a LNV effective lagrangian of the form 

£(M3/3 = % °* 0* = ^ Tiq ^ r 2qeT 3 e c , (18) 

i 

where we have suppressed the flavor and Dirac structures (a complete list of the dimension nine operators 
Oi can be found in Ref. |64j). 

For the CLFV interactions, one has operators of dimension six, 

£clfv = £^0,, (19) 

i 

The CLFV operators relevant to our analysis are of the following type (along with their analogues with 
L^R): 

O aL = j^I~ L a, u ipi jL + h.c. , O eL = I~ L l) L if L £ mL , O iq = ^ qT q q . (20) 
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Operators of the type O a are typically generated at one-loop level, hence our choice to explicitly 
display the loop factor 1/(Att) 2 . On the other hand, in a large class of models, operators of the type Oi 
or Oe q are generated by tree level exchange of heavy degrees of freedom. With the above choices, all 
non-zero Cj are nominally of the same size, typically the product of two Yukawa-like couplings or gauge 
couplings (times flavor mixing matrices). These operators are analogous to the operators in Eqs. ^ 
and ^ but normalized differently in order to follow the notation of Ref . [66] . 

Whenever the operators 0(l,r and/or Og q appear at tree-level in the effective theory, they lead to an 

enhancement factor log ^^r^ in the LFV processes involving charged muons [i — \ e (muon conversion) 

or fi — y eee. That logarithmic enhancement is absent in the \i — > process. 

Based on our analysis in Ref. [66J we can formulate the following approximate rules ( a diagnostic 
tool) regarding the relative size of the branching ratios for fi — > e and \i — > ej and the relation of that 
size to the mechanism of the Ov/3/3 decay: 

1. Observation of both the CLFV muon processes // — > e and \i — > with relative branching ratios 
~ 10~ 2 implies, under generic conditions, that ~ (m^) 2 . Hence the relation of the 0^(3(3 
lifetime to the absolute neutrino mass scale is straightforward. 

2. On the other hand, observation of CLFV muon processes with relative branching ratios 3> 1CT 2 
(presumably caused by the log enhancement factor ) could signal non-trivial LNV dynamics 
at the TeV scale, whose effect on 0^(3(3 has to be analyzed on a case by case basis. Therefore, in 
this scenario no definite conclusion can be drawn based on LFV rates. 

3. Non-observation of CLFV in muon processes in forthcoming substantially more sensitive experi- 
ments would imply either that the scale of non-trivial CLFV and LNV is above a few TeV, and 
thus T^pp ~ (m^) 2 , or that any TeV-scale LNV is approximately flavor diagonal. 

The above statements were illustrated using two explicit cases: the Minimal Supersymmetric Stan- 
dard Model with R-parity violation (RPV-SUSY) and the Left-Right Symmetric Model (LRSM) in Ref. 



It is likely that the basic mechanism at work in these illustrative cases is generic: low scale LNV 
interactions (AL = ±1 and/or AL = ±2), which in general contribute to Oz//3/3, also generate sizable 
contributions to fi — > e conversion, thus enhancing this process over fi — > ej (see also discussion following 



Eq. (10) in Section IID. 



3.6 0v(5(5 nuclear matrix elements 

Double beta decay in both 2v and 0u modes can exist because even-even nuclei are more bound than 
the neighboring odd-odd nuclei. This extra binding is a consequence of pairing between like nucleons. 
In nonmagic systems neutrons and/or protons form + pairs and the corresponding Fermi level becomes 
diffuse over the region with the characteristic size ~ pairing gap A. This opens more possibilities for 
nn — > pp transitions. The calculated matrix elements M 0v increase when the gap A increases. The 
(unrealistic in real nuclei) situation of pure paired nuclear system (only seniority states) would have 
very large M 0v . 

However, in real nuclei opposite tendencies are also present. Real nuclei have admixtures of the 
"broken pair" states, or in the shell model language, states with higher seniority. These states are 
present because other parts of the nucleon interaction exist, in particular the neutron-proton force. It 
is illustrative to characterize such states by the angular momentum J of the neutron pair that is in the 
(3(3 decay transformed into the proton pair with the same J . While the pairing parts J = are large 
and positive, the J ^ parts are negative and their sum is essentially as large as the J = piece. The 
severe cancellation between these two tendencies, as a consequence of the corresponding components of 
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the residual interaction, is present in most methods of evaluating the M 0u . That makes their accurate 
determination challenging. 

The Qv/3/3 decay operators depend on the internucleon distance ri 2 due to the neutrino potential 
H(r,E) (which is essentially the Fourier transform of the light neutrino propagator). Obviously, the 
range of r 12 is restricted from above by r 12 < 2R. From the form of H(r) ~ R/r one could, naively, 
expect that the characteristic value of r 12 is the typical distance between nucleons in the nucleus, namely 
fi2 ~ R. However, that is not true, in reality only much smaller values of ri 2 < 2-3 fm or equivalently 
larger values of the momentum transfer q are relevant. The qualitatively same conclusion was reached 
in both widely used approximate methods, the Quasiparticle Random Phase Approximation (QRPA) 
[67] and in the Nuclear Shell Model (NSM) [68]. Examples of the functions C 0u (r) are shown in Fig. [7] 
for three representative nuclei. As the lower panel demonstrates, the cancellation between the "pairing" 
{J = 0) and "broken pairs" {J ^ 0) is essentially complete for r i2 > 2-3 fm 

To see how that conclusion is obtained, we define the function C{r) 

C°J T (r) = (f\^ikffi-^kr^8(r-r lk )H(r lk ,E)\i) , (21) 
Obviously, this function is normalized by 

POO 

M% T = / C% T {r)dr . (22) 
Jo 

The definition above is valid for the Gamow- Teller part of M 0v ', which is, however, its most important 
part. In Fig. [7] the full C 0u (r) is plotted (analogous figure appeared in [67]). Note that due to the 
relatively short range of C 0u (r) proper treatment of the nucleon short-range repulsion, nucleon finite 
size, as well as of the recoil order parts of the nucleon weak current (induced pseudoscalar, weak 
magnetism) are necessary. 

Variety of methods has been used for evaluation of the OvftP nuclear matrix elements for the simplest 
scenario where the decay rate is proportional < mpp > 2 . They differ in their choice of the valence space, 
interaction hamiltonian and the ways the corresponding equations of motion are solved. Let us stress 
that an exact u ab initio" , i.e. without approximations, calculation of M 0v for the candidate nuclei is 
impossible at the present time. 

The nuclear shell model (NSM) is, in principle, the method that seems to be well suited for this task. 
In it, the valence space consists of just few single particle states near the Fermi level. With interaction 
that is based on the realistic nucleon-nucleon force, but renormalized slightly to describe better masses, 
energies and transition probabilities in real nuclei, all possible configurations of the valence nucleons 
are included in the calculation. The resulting states have not only the correct number of protons and 
neutrons, but also all relevant quantum numbers (angular momentum, isospin, etc). 

The quasiparticle random phase approximation (QRPA) and its renormalized version (RQRPA) is 
another method often used in the evaluation of M 0u . In it, the valence space is not restricted and 
contains at least two full oscillator shells, often more than that. On the other hand, only selected 
simple configurations of the valence nucleons are used. The basis states have broken symmetries in 
which particle numbers, isospin, and possibly angular momentum are not good quantum numbers but 
conserved only on average. After the equations of motion are solved, some of the symmetries are 
partially restored. The RQRPA partially restores the Pauli principle violation in the resulting states. 
Thus, in certain sense the NSM and QRPA are complementary methods. 

The IBM-2 method uses the microscopic interacting boson model to evaluate M 0u [6H] ■ In IBM-2 one 
begins with correlated S (angular momentum 0) and D (angular momentum 2) pairs of identical nucleons 
and includes the effect of deformation through the bosonic neutron-proton quadrupole interaction. The 
method describes well the low lying states, the electromagnetic transitions between them and the 
two- nucleon transition rates in spherical and strongly deformed nuclei. Even though the method was 



19 




Figure 7: The dependence on r 12 for 76 Ge, 100 Mo and 130 Te evaluated in QRPA. The upper 
panel shows the full matrix element, and the lower panel shows separately the "pairing" {J 
= 0) and "broken pair" {J ^ 0) contributions. The integrated M 0u are 5.4, 4.5, and 4.1 for 
the three nuclei. 



originally considered as an approximation of the nuclear shell model, the resulting M 0u are, rather 
surprisingly, close to the RQRPA results and noticeably larger than the NSM ones. At the present time 
it is not possible to evaluate the M 2v matrix elements within the IBM-2 method, except in the closure 
approximation. 

In Ref. [70] the generating coordinate method (GCM or EDF) was employed. Large single particle 
space was used (11 shells) with the well established Gogny D1S energy density functional. The initial and 
final many-body wave functions are represented as combinations of the particle number N, Z projected, 
I = + axially symmetric states with different intrinsic deformations. Again, like in IBM-2 method, 
the treatment of the odd-odd nuclei, and thus also of the M 2u matrix elements is impossible in GCM. 

Finally, the PHFB method, used in Ref. [71], uses the projected Hartree-Fock-Bogolyubov wave 
function. 

The results of all of these five methods for the most important candidate nuclei are shown in Fig. [8] 
(NSM is from [HE], IMB-2 from [69], EDF from [70], RQRPA from [72] and [73]) and PHFB from [71]. 
There one can see that, common to all displayed methods, the predicted M 0v nuclear matrix elements 
vary relatively smoothly, with the mass number A, unlike the experimentally determined M 2u matrix 
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Figure 8: Dimensionless Ov/3/3 nuclear matrix elements for selected nuclei evaluated using a 
variety of indicated methods. For references see text. 



elements displayed in Fig. [6] that vary strongly with A. The RQRPA, IBM-2 and GCM methods are 
in a crude agreement with each other, and predict slow decrease of M 0u with increasing A. The M 0u 
evaluated in NSM are essentially constant with A and noticeably smaller than those from the other 
methods, particularly in the lighter nuclei. 

As stressed earlier, when heavy particles of any kind mediate the Ov/3/3 decay, we are dealing with a 
six fermion vertex, representing extremely short range operator. Traditionally, the short range nucleon- 
nucleon repulsion, which would heavily suppress the matrix elements in that case, is overcome by 
introducing the dipole type nucleon form factor (see [74j). The neutrino potential will be then of the 
form 



heavyu \ 
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Ml 



(27T) 



Mi 



M| + <f 



M A R 
48 ( 



-M A r 
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(23) 



Such potential will have the range 1/M A ~ 1/GeV and will be much less affected by the short range 
correlations. The disadvantage of the form factor modeling is that the error introduced by such approx- 
imation is very difficult to estimate. 

However, as pointed out in Ref. [75] and analyzed in detail in [75] using the Effective Field Theory 
(EFT) approach, the more effective and fundamentally leading effect is to replace the NNNNee six 
fermion vertex by the nnee vertex. Concrete application of the pion exchange mechanism suggests the 
dominance of the irnee over the short range nucleon only vertex by a factor of 10 - 30 at the level of 
the nuclear matrix element. This remains, to some extent, so far largely unexplored issue. 



21 



3.7 Possible existence of light sterile neutrinos 



Most particle physics models of neutrino mass contain additional, so far unobserved, gauge singlet 
neutral fermions, commonly identified as sterile or right-handed neutrinos. Such Vr could be very 
heavy, such as in the orthodox Type I seesaw, and hence unobservable. In other models (e.g. in seesaw 
Type II and III) such sterile neutrinos have ~TeV masses, and their mixing with the three active 
neutrinos is again expected to be very small. Existence of sterile neutrinos in models of this type is 
needed and not surprising; without them proper explanation of the masses of the observed light active 
neutrinos would be impossible. 

Remarkably, recently there appeared several empirical indications, so far of somewhat limited sta- 
tistical significance and mutual consistency, that very light (of mass 0(leV)) sterile neutrinos might 
exist, that moreover mix noticeably with the flavor neutrinos. Since these hints appear in a variety 
of independent experiments, it is worthwhile to analyze their possible effect on the problem discussed 
here, namely on the OvfiP decay. One has to keep in mind, however, that if the existence of such sterile 
neutrinos is confirmed, it would be necessary to reformulate the underlying models and find convincing 
physics reasons for their existence and for the magnitude of their mixing. 

The well known LSND observation of the v e appearance interpreted as the — > v e oscillation [77] 
and the results of the MiniBooNE experiment at similar L/E v \J3\ are consistent with oscillation at 
a Am 2 ~ 1 eV 2 . The so-called reactor anomaly, a deficit of the measured reactor i> e flux at L < 100 
m compared to the reevaluated flux prediction [79"| |8*U] . as well as the analysis of the calibration of 
Gallex-SAGE solar neutrino detectors with the radioactive sources [81 J have significance of ~ 3a and 
point toward a similar interpretation. On the other hand, a number of short baseline oscillation searches 
have yielded null results, and a consistent and convincing picture of all of these effects has not been 
achieved as yet. 

Possible existence of ~ 1 eV sterile neutrino on the effective neutrino Majorana mass < mpp > 
explored in the 0v(3(3 decay would be profound. < mpp > depends on the mass eigenstates neutrinos Vi 
with masses rrii as 

< mpp >= |S i |f/ ei | 2 m i exp(zai)| , (24) 

where U e i are the entries of the first row of the neutrino mixing matrix and are the so-called Majorana 
phases. 



When the summation in Eq. (24) goes only over the three known light Vi, as it is well known, in 
the case of the inverted hierarchy there is a lower limit < mpp >> 20 meV. If the inverted hierarchy 
is shown to be the correct pattern of neutrino masses experimentally (not by 0v(3(3 decay), the above 
statement would imply that reaching the sensitivity to < mpp >~ 20 meV would represent a definitive 
test of the Majorana nature of neutrinos. The situation changes drastically when the existence of a 
fourth, necessarily sterile, neutrino is assumed. The lower limit of < mpp > disappears in that case; 
the whole region up to < mpp > ~ 0.1 eV is now allowed. The best fit region for the normal hierarchy 
is changed substantially as well. 



3.8 Where are we now 

As pointed out in connection with Fig. [5j there has been a decade long hiatus in the improvements of 
the sensitivity of the searches for 0u/3/3 decay when expressed as the upper limits of < mpp >. The 
2001 result of the Heidelberg-Moscow experiment, 1.9 x 10 25 years of the half-life for 76 Ge, corresponds 
to < mpp > in the range 0.25-0.5 eV depending on the calculated value of the corresponding nuclear 
matrix element. Later, a subset of that collaboration reanalyzed the experimental result and claimed 
[82] a positive observation of the 0v(3(3 decay with Ti/ 2 = 2.23^g;3g x 10 25 years, implying < mpp > = 
0.32 ± 0.03 eV using their preferred matrix elements [86J. That result remains unconfirmed. 
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Figure 9: Relation between the Ou/3/3 half-lives of Xe and Ge. The claim [52] is rep- 
resented by the grey horizontal band along with the best limit for 76 Ge [83] . The 136 Xe 
limits from [84J and [85] are shown by the vertical lines. Different theoretical calculations 
are represented by the diagonal lines, (GCM from [70], NSM from [68], IBM-2 from [69] . 
RQRPA-1 from [72], and QRPA-2 from [SB]) with tick marks indicating the corresponding 
< mpp > values. 



Very recently two experimental searches for the Op/3/3 decay of 136 Xe published their results. In Fig. [9] 
these recent results are compared with those for 76 Ge. At 90% CL the EXO-200 result corresponds to 
< mpp > of less than 0.14 - 0.38 eV. At 68% (90%) CL it contradicts the claim [52] for all (most) matrix 
element element calculations (analogous figure appears in [85]) . 

Both competing 136 Xe experiments [541 155] are continuing their runs and substantial improvement 
of their half-life sensitivities both through accumulating statistics and improving performance of their 
detectors is expected. Moreover, other experiments of comparable size and sensitivity are expected 
to begin their operation soon. The GERDA experiment [ST] is using 76 Ge again, and will definitely 
confirm or reject the claim [82]. The CUORE experiment [55] is using 130 Te and will complement the 
results obtained with 76 Ge and 136 Xe. There is, therefore, a realistic chance that the upper edge of the 
so-called degenerate mass region ((m^g) ~ 0.05 - 0.1 eV) will be reached reasonably soon. 

Even more ambitious proposals for ton (or multi-ton) size experiments to search for the neutrinoless 
(3/3 decay exist. If the projected performance (sensitivity, background suppression) could be achieved, 
and the corresponding funding secured, we should, perhaps within a decade, reach the boundary of 
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( m /3^) ~ 20 meV. It is very important, essentially imperative, to continue the search for the 0^/3/3 
decay, by a variety of techniques and using several candidate nuclei as sources. Only this way, if the 
OvfiP decay is actually observed, we will have do doubts of its discovery. 

To extend the sensitivity beyond {mpp) ~ 20 meV, the region where only the normal hierarchy 
is present, appears to be extremely challenging. No concrete plans exist for such experiments at the 
present time. 

Above, we used the effective Majorana mass < mpp > to indicate the sensitivity of the various 
searches for the Ov/3/3 decay. However, as pointed out already, this is not the only possible mecha- 
nism. Decay mediated by various heavy particles is also a possibility and we have to keep it in mind 
when comparing the 0^/3/3 decay neutrino mass limits with the limits obtained in the direct kinematic 
experiments and from cosmology/astrophysics. 

As we tried to stress repeatedly in this article, the fundamental importance of the Qv/3/3 decay is the 
test of the total lepton number conservation. If discovered, it would be a crucial step in the possible 
formulation of the 'New Standard Model'. 
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